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^ ' Abstract 
. 

I calculate the classical effects induced by an isotropic mass loss 
■ M/M of a body on the orbital motion of a test particle around it; 

the present analysis is also valid for a variation G/G oi the Newto- 
nian constant of gravitation. I perturbatively obtain negative secular 
rates for the osculating semimajor axis a, the eccentricity e and the 
mean anomaly A4, while the argument of pericenter w does not un- 
^T) [ dergo secular precession; the longitude of the ascending node il and 

the inclination i remain unchanged as well. The anomalistic period is 
different from the Keplerian one, being larger than it. The true orbit, 
- instead, expands, as shown by a numerical integration of the equations 

I of motion in Cartesian coordinates; in fact, this is in agreement with 

the seemingly counter-intuitive decreasing of a and e because they only 
refer to the osculating Keplerian ellipses which approximate the tra- 
jectory at each instant. By assuming for the Sun M/M = —9 x 10^^"* 
^j^l yr~^ it turns out that the Earth's perihelion position is displaced out- 

ward by 1.3 cm along the fixed line of apsides after each revolution. 
By applying my results to the phase in which the radius of the Sun, 
. already moved to the Red Giant Branch of the Hertzsprung-Russell 

5-H ' Diagram, will become as large as 1.20 AU in about 1 Myr, I find that 

^ the Earth's perihelion position on the fixed line of the apsides will 

increase by « 0.22 - 0.25 AU (for M/M = -2 x 10"'^ yr'^); other 
researchers point towards an increase of 0.37 — 0.63 AU. Mercury will 
be destroyed already at the end of the Main Sequence, while Venus 
should be engulfed in the initial phase of the Red Giant Branch phase; 
the orbits of the outer planets will increase by 1.2 — 7.5 AU. Simultane- 
ous long-term numerical integrations of the equations of motion of all 
the major bodies of the solar system, with the inclusion of a mass-loss 
term in the dynamical force models as well, are required to check if 
the mutual N-body interactions may substantially change the picture 
analytically outlined here, especially in the Red Giant Branch phase 
in which Mercury and Venus may be removed from the integration. 

Keywords: gravitation, stars: mass-loss, celestial mechanics 
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1 Introduction 



I deal with the topic of determining the classical orbital effects induced by 
an isotropic variation M /M of the mass of a central body on the motion 
of a test particle; my analysis is also valid for a change G/G of the Newto- 
nian constant of gravitation. This problem, although interesting in itself, is 
not only an academic one because of the relevance that it may have on the 
ultimate destiny of planetary companions in many stellar systems in which 
the host star experiences a mass loss, like our Sun [T]. With respect to 
this aspect, my analysis may be helpful in driving future researches towards 
the implementation of long-term N-body simulations including the tempo- 
ral change of GM as well, especially over timescales covering paleoclimate 
changes, up to the Red Giant Branch (RGB) phase in which some of the 
inner planets should be engulfed by the expanding Sun. Another problem, 
linked to the one investigated here, which has recently received attention is 
the observationally determined secular variation of the Astronomical Unit 
[21 [3l m [5] . Moreover, increasing accuracy in astrometry pointing towards 
microarcsecond level [6], and long-term stability in clocks [7] require to con- 
sider the possibility that smaller and subtler perturbations will be soon de- 
tectable in the solar system. Also future planetary ephemerides should take 
into account M/M. Other phenomena which may show connections with 
the problem treated here are the secular decrease of the semimajor axes of 
the LAGEOS satellites, amounting to 1.1 mm d-\ [8] and the increase of 
the lunar orbit's eccentricity [9|. However, a detailed analysis of all such 
issues is beyond the scope of this paper. 

Many treatments of the mass loss-driven orbital dynamics in the frame- 
work of the Newtonian mechanics, based on different approaches and laws 
of variation of the central body's mass, can be found in literature; see, e.g., 
[inillllliaiiainillllllSlinillHlElllland references therein. 

The plan of the paper is as follows. Section [2] is devoted to a theoretical 
description of the phenomenon in a two-body scenario. By working in the 
Newtonian framework, I will analytically work out the changes after one 
orbital revolution experienced by all the Keplerian orbital elements of a test 
particle moving in the gravitational field of a central mass experiencing a 
variation of its GM linear in time. Then, I will clarify the meaning of the 
results obtained by performing a numerical integration of the equations of 
motion in order to visualize the true trajectory followed by the planet. Con- 
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cerning the method adopted, I will use the Gauss perturbation equations 
|19l 120] . which are valid for generic disturbing accelerations depending on 
position, velocity and time, the "standard" Keplerian orbital elements (the 
Type I according to, e.g., [H]) with the eccentric anomaly E as "fast" angu- 
lar variable. Other approaches and angular variables like, e.g. the Lagrange 
perturbation equations [191 [2Q], the Type II orbital elements IT6| and the 
mean anomaly A4 could be used, but, in my opinion, at a price of major 
conceptual and computational difficultie^. With respect to possible connec- 
tions with realistic situations, it should be noted that, after all, the Type I 
orbital elements are usually determined or improved in standard data reduc- 
tion analyses of the motion of planets and (natural and artificial) satellites. 
Instead, my approach should, hopefully, appear more transparent and easy 
to interpret, although, at first sight, some counter- intuitive results concern- 
ing the semimajor axis and the eccentricity will be obtained; moreover, for 
the chosen time variation of the mass of the primary, no approximations 
are used in the calculations which are quite straightforward. However, it is 
important to stress that such allegedly puzzling features are only seemingly 
paradoxical because they will turn out to be in agreement with numerical 
integrations of the equations of motion, as explicitly shown by the Figures 
depicted. Anyway, the interested reader is advised to look also at [Mj for 
a different approach. In Section [3] I will apply my results to the future 
Sun-Earth scenario and to the other planets of the solar system. Section |4] 
summarizes my results. 



2 Analytical calculation of the orbital effects by 

By defining 



/i = GM 



at a given epoch to, the acceleration of a test particle orbiting a central body 
experiencing a variation of ^ is, to first order in t — to, 



ii 



1 + 



{t - to) 



(2) 



with fi = fi\t=to- A "^ill be assumed constant throughout the temporal 
interval of interest At = t — to, as it is, e.g., the case for most of the 



^ Think, e.g., about the cumbersome expansions in terms of the mean anomaly and the 
Hansen coefficients, the subtleties concerning the choice of the independent variable in the 
Lagrange equations for the semimajor axis and the eccentricity |19| . 
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remaining lifetime of the Sun as a Main Sequence (MS) star [I]. Note that fi 
can, in principle, be due to a variation of both the Newtonian gravitational 
constant G and the mass M of the central body, so that 

+ ^ (3) 

Moreover, while the orbital angular momentum is conserved, this does not 
happen for the energy. 

By limiting ourselves to realistic astronomical scenarios like our solar 
system, it is quite realistic to assume that 



{t - to) < 1 (4) 



over most of its remaining lifetime: indeed, since M/M is of the order o^l 
10"^'^ yr~^ for the Sun the condition of eq. is satisfied for the 
remaining ~ 7.58 Gyr before the Sun will approach the RGB tip in the 
Hertzsprung- Russell Diagram (HRD). Thus, I can treat it perturbatively 
with the standard methods of celestial mechanics. 

The unperturbed Keplerian ellipse at epoch t^, assumed coinciding with 
the time of the passage at perihelion tp, is characterized by 

r = a(l — ecosE'), 
dt = {^=^f^) dE, 



cos/ 



cos hi—e 
1— ecos E ' 



1— ecosi? ' 



sin/ 

where a and e are the semimajor axis and the eccentricity, respectively, which 
fix the size and the shape of the unchanging Keplerian orbit, n = \J ^/a^ 
is its unperturbed Keplerian mean motion, / is the true anomaly, reckoned 
from the pericentre, and E is the eccentric anomaly. Eq. ([5|) characterizes 
the path followed by the particle for any t > tp ii the mass loss would 
suddenly cease at tp. Instead, the true path will be, in general, different from 



About 80% of such a mass-loss is due to the core nuclear burning, while the remaining 
20% is due to average solar wind. 

^The age of the present-day MS Sun is 4.58 Gyr, counted from its zero-age MS star 
model pp. 
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a closed ellipse because of the perturbation induced by /i and the orbital 
parameters of the osculating ellipses approximating the real trajectory at 
each instant of time will slowly change in time. 



2.1 The semimajor axis and the eccentricity 

The Gauss equation for the variation of the semimajor axis a is [19^ [20] 



da 
lit 



n 



eAr sin f + A^- (- 



(6) 



where Ar and A^ are the radial and transverse, i.e. orthogonal to the di- 
rection of f , components, respectively, of the disturbing acceleration, and 
p = a(l — e^) is the semilatus rectum. In the present case 



(t-tp), 



(7) 



i.e. there is an entirely radial perturbing acceleration. For /i < 0, i.e. a 
decrease in the body's GM, the total gravitational attraction felt by the 
test particle, given by eq. is reduced with respect to the epoch tp. In 
order to have the rate of the semimajor axis averaged over one (Keplerian) 
orbital revolution eq. ^ must be inserted into eq. ([6]), evaluated onto 
the unperturbed Keplerian ellipse with eq. ([5]) and finally integrated over 
ndt/2-K from to 2-k because n/27r = l/P^^^P (see below). Note that, from 
eq. ([5]), it can be obtained 



E — e sin E 



n 



As a result, I hav^ 



da 



- - a 



2tt 



[E — e sin E) sin E 
(1 — ecosE')^ 



dE 



1 



^^a. 



(8) 



(9) 



Note that if ^ decreases a gets reduced as well: (d) < 0. This may be 
seemingly bizarre and counter-intuitive, but, as it will be shown later, it is 
not in contrast with the true orbital motion. 

The Gauss equation for the variation of the eccentricity is [191 EQ] 



de 



na 



Ay. sin / + At 



1 / r 
cos / + - 1 - - 
e V a 



(10) 



^Recall that the integration is taken over the unperturbed Keplerian elhpse: that is 
why a and e are kept out of the integral in eq. ([9]) and in the following averages. 
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For A = Ar, it reduces to 



de f 1 — da 



dt \ 2ae J dt 
so that 

also the eccentricity gets smaller for /i < 0. 

As a consequence of the found variations of the osculating semimajor 
axis and the eccentricity, the osculating orbital angular momentum per unit 
mass, defined by = fia{l — e^), remains constant: indeed, by using eq. 
([9]) and eq. ([T2|), it turns out 

dL"^ \ 

/i (d) (1 - e^) - 2/iae (e) = 0. (13) 



dt 

The osculating total energy £ = —fi/2a decreases according to 

Moreover, the osculating Keplerian period 

pKcp ^ 2Tr^, (15) 

which, by definition, yields the time elapsed between two consecutive per- 
ihelion crossings in absence of perturbation, i.e. it is the time required to 
describe a fixed osculating Keplerian ellipse, decreases according to 

/dP^\ ^ 3 Kep(«) ^ 6^eA faV/' 
\ dt / 2 a (1 - e) V/V ■ 

As I will show, also such a result is not in contrast with the genuine orbital 
evolution. 

2.2 The pericentre, the node and the inchnation 

The Gauss equation for the variation of the pericentre to is |19l [20] 



du} \/l — e 



2 



dt nae 



-ArCosf + Ar{l + ^] sin/ 



COSZ-, (17) 
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where i and il. are the the inchnation and the longitude of the ascending 
node, respectively, which fix the orientation of the osculating ellipse in the 
inertial space. Since (Kl/dt and di/dt depend on the normal component 
Ajj of the disturbing acceleration, which is absent in the present case, and 
A = Ar, I have 

dt I 27re \li)h (l-ecos£;)2 ^ ' 

the osculating ellipse does not change its orientation in the orbital plane, 
which, incidentally, remains fixed in the inertial space because Ay = and, 
thus, dVL/dt = di/dt = 0. 



2.3 The mean anomaly 

The Gauss equation for the mean anomaly A4, defined as = n(t — tp), 
[El [20] is 

dM 2 r r -(duj .dn\ , , 

— =n A^_-v^l-e2 — +cos^— . (19) 

dt na a \ at at J 

It turns out that, since 

2 r 2/1 

Ar-dt = -^{E - e sin E)dE, (20) 

na a n'^a'^ 

then 

^ n + 27r^ ; 21 



dt / \ /i ^ 

the mean anomaly changes uniformly in time at a slower rate with respect 
to the unperturbed Keplerian case for /i < 0. 



2.4 Numerical integration of the equations of motion and 
explanation of the seeming contradiction with the ana- 
lytical results 

At first sight, the results obtained here may be rather confusing: if the 
gravitational attraction of the Sun reduces in time because of its mass loss 
the orbits of the planets should expand (see the trajectory plotted in Figure 
m numerically integrated with MATHEMATICA), while I obtained that the 
semimajor axis and the eccentricity undergo secular decrements. Moreover, 
I found that the Keplerian period P^'^p decreases, while one would expect 
that the orbital period increases. 
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y, AU 




Figure 1: Black continuous line: true trajectory obtained by numerically 
integrating with MATHEMATICA the perturbed equations of motion in 
Cartesian coordinates over 2 yr; the disturbing acceleration of eq. ([2]) has 
been adopted. The planet starts from the perihelion on the x axis. Just for 
illustrative purposes, a mass loss rate of the order of 10~^ y^~^ has been 
adopted for the Sun; for the planet initial conditions corresponding to a = 1 
AU, e = 0.8 have been chosen. Red dashed line: unperturbed Kepler ian 
ellipse at t = Iq = tp. Blue dash-dotted line: osculating Keplerian ellipse 
after the first perihelion passage. As can be noted, its semimajor axis and 
eccentricity are clearly smaller than those of the initial unperturbed ellipse. 
Note also that after 2 yr the planet has not yet reached the perihelion as 
it would have done in absence of mass loss, i.e. the true orbital period is 
longer than the Keplerian one of the osculating red ellipse. 
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In fact, there is no contradiction, and my analytical results do yield us 
realistic information on the true evolution of the planetary motion. Indeed, 
o, e and P^^p refer to the osculating Keplerian ellipses which, at any instant, 
approximate the true trajectory; it, instead, is not an ellipse, not being 
bounded. Let us start at tp from the osculating pericentre of the Keplerian 
ellipse corresponding to chosen initial conditions: let us use a heliocentric 
frame with the x axis oriented along the osculating pericentre. After a 
true revolution, i.e. when the true radius vector of the planet has swept 
an angular interval of 27r, the planet finds itself again on the X axis, but 
at a larger distance from the starting point because of the orbit expansion 
induced by the Sun's mass loss. It is not difficult to understand that the 
osculating Keplerian ellipse approximating the trajectory at this perihelion 
passage is oriented as before because there is no variation of the (osculating) 
argument of pericentre, but has smaller semimajor axis and eccentricity. 
And so on, revolution after revolution, until the perturbation theory can be 
applied, i.e. until — tp) « 1. In Figured] the situation described so 

far is qualitatively illustrated. Just for illustrative purposes I enhanced the 
overall effect by assuming ~ 10"^ yr"^ for the Sun; the initial conditions 
for the planet correspond to an unperturbed Keplerian ellipse with a = 1 AU, 
e = 0.8 with the present-day value of the Sun's mass in one of its foci. It is 
apparent that the initial osculating red dashed ellipse has larger a and e with 
respect to the second osculating blue dash-dotted ellipse. Note also that the 
true orbital period, intended as the time elapsed between two consecutive 
crossings of the perihelion, is larger than the unperturbed Keplerian one of 
the initial red dashed osculating ellipse, which would amount to 1 yr for the 
Earth: indeed, after 2 yr the planet has not yet reached the perihelion for 
its second passage. 

Now, if I compute the radial change Ar{E) in the osculating radius vector 
as a function of the eccentric anomaly E I can gain useful insights concern- 
ing how much the true path has expanded after two consecutive perihelion 
passages. Prom the Keplerian expression of the Sun-planet distance 



one gets the radial component of the orbital perturbation expressed in terms 
of the eccentric anomaly E 



r = a(l — e cos E) 




Ar{E) = {1 - e cos E) Aa - a cos E Ae + ae sin E AE; 



(23) 
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it agrees with the results obtained in, e.g., [21]. Since 

zae / A* 1 / sin cos _E 



Ae 



2ae I /£ 
n \fi 



1— e cos i? 



sinE—EcosE 
1— e cos E 



(24) 



with 



C{E) 



_E^+2e(cos g-1) 
1— ecosi? ' 

1-e 



l+e-(l+e) cos_E-Esin_E 
(l-ecos£;)2 



(1— e^) sini?(sin £— ecos E) 



(25) 



(1— e cos E)'-' 



it follows 



Ar(i^)= « (^)[p(ii;) + ^(ii;)], 



(26) 



with 



P(S) 



HE) 



-2(sin S - cos S) + _ (l-e^)sin^i;(sini; ecosi.) - 

/ 1— ecosi? (1— ecosi?)^ 



1— e cos i? 



COS £'(sin E — E cos -E) + sin E 



l+e- (l+e) cos E-E sin E 
1—e cos E 



(27) 

It turns out from eq. (|26p and eq. (j27p that, for > 0, Ar{E) never 
vanishes; after one orbital revolution, i.e. after that an angular interval of 
2tt has been swept by the (osculating) radius vector, a net increase of the 
radial (osculating) distance occurs according tcH 



Ar(27r) - Ar(0) = Ar(27r) = -—a ( - ) (1 - e). 



(28) 



This analytical result is qualitatively confirmed by the differenc^ Ar{t) 
between the radial distances obtained from the solutions of two numerical 
integrations of the equations of motion over 3 yr with and without /t///; the 
initial conditions are the same. For illustrative purposes I used o = 1 AU, 
e = 0.01, fi/fi = —0.1 yr~^. The result is depicted in Figure [2j Note also 



■■^According to eq. ^ and eq. Ar(0) = 0. 

^Strictly speaking, Ar and the quantity plotted in Figure [2] are different objects, but, 
as the following discussion will clarify, I can assume that, in practice, they are the same. 
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Ar (AU) 




Figure 2: Difference Ar{t) between the radial distances obtained from the 
sohitions of two numerical integrations with MATHEMATICA of the equa- 
tions of motion over 3 yr with and without fi/n; the initial conditions are 
the same. Just for illustrative purposes a mass loss rate of the order of 
— 10~^ yr~^ has been adopted for the Sun; for the planet initial conditions 
corresponding to a = 1 AU, e = 0.01 have been chosen. The cumulative 
increase of the Sun-planet distance induced by the mass loss is apparent. 
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True orbit 



Unperturbed ellipse 
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\ 
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Line of the 
apsides 



Figure 3: Radial and transverse perturbations Ar and At of the Keplerian 
radius vector (in blue); the presence of the transverse perturbation At 
makes the real orbit (in red) lagging behind the Keplerian one. 

that eq. (p6|) and eq. ([27|) tell us that the shift at the aphelion is 



in agreement with Figure [T] where it is 4.5 times larger than the shift at the 
perihelion. 

Since Figure [1] tells us that the orbital period gets larger than the Kep- 
lerian one, it means that the true orbit must somehow remain behind with 
respect to the Keplerian one. Thus, a negative perturbation Ar in the 
transverse direction must occur as well; see Figure [3l 

Let us now analytically compute it. According to [21j, it can be used 




(29) 



Ar 



asinE 



+ a^/l-e^ AE + r{Auj + Ail cos i). 



(30) 
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By recalling that, in the present case, AO, = and using 

■ 1 + e - (1 + e) cos E' - sin 



Auj 



ne 



1 — e cos E 



(31) 



it is possible to obtain from eq. (j24p and eq. ([3T 



n V/uy (1 — ecos E) 



[GiE) + niE) + 1{E) + JiE) + }CiE)] 



(32) 



with 



g{E) 
n{E) 

X{E) 
J{E) 
)C{E) 



sin E{E cos — sin E) , 

(i-ecosE) [(i + e)(cosE- l) + EsinE], 

£;2 + 2e(cosE-l), 

"(l-e^)(ecos£;-sin£;)" 



(33) 



sin£' 



1— e cos i? 



(l+e)(l-ecos E)-EsinE 
1— e cos £ 



It turns out from eq. ([32]) and eq. (j33|) that, for E > 0, At{E) never 
vanishes; at the time of perihelion passage 



Ar(27r) - Ar(0) = a ^ 



4tt^ f fi\ l + e 



n \fi J \ 1 — e 



< 0. 



(34) 



This means that when the Keplerian path has reached the perihelion, the 
perturbed orbit is still behind it. Such features are qualitatively confirmed 
by Figured) 

From a vectorial point of view, the radial and transverse perturbations 
to the Keplerian radius vector r yield a correction 



A = Ar f + Ar f. 



so that 



The length of A is 



pert 



r + A. 



A{E) = v/Ar(E)2 + Ar(E)2; 



(35) 
(36) 

(37) 
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eq. ([28]) and eq. ([32]) tell us that at perihelion it amounts to 



A(2^) = Ar(2^)^/l + 47r2i^-l^. (38) 



(1 + e 
(1-e) 

The angle ^ between A and r is given by 

*^^^(^) = A^= ^^^^ 

at perihelion it is 

tan^(2vr) = -2vr^^^±|^, (40) 

i.e. ^ is close to —90 deg; for the Earth it is —81.1 deg. Thus, the difference 
5 between the lengths of the perturbed radius vector rpert and the Keplerian 
one r at a given instant amounts to about 

5«AcosC; (41) 

in fact, this is precisely the quantity determined over 3 yr by the numerical 
integration of Figure [2j At the perihelion I have 



5 = Ar (2^) W 1 + 47r2 l\ ^ f „ cos g; (42) 



since for the Earth 



Wl + 47r2ii-t4l cos ^ = 1.0037, (43) 
V (1-e)^ 

it holds 

5 ^ Ar(27r). (44) 

This explains why Figure [2] gives us just Ar. 

Concerning the observationally determined increase of the Astronomical 
Unit, more recent estimates from processing of huge planetary data sets by 
Pitjeva point towards a rate of the order of 10"^ m yr~^ |22t I23j . It may 
be noted that my result for the secular variation of the terrestrial radial 
position on the line of the apsides would agree with such a figure by either 
assuming a mass loss by the Sun of just —9 x 10~^^ yr"^ or a decrease of the 
Newtonian gravitational constant G/G ~ — 1 x 10~^^ yr^^- Such a value for 
the temporal variation of G is in agreement with recent upper limits from 
Lunar Laser Ranging [24J G/G = (2 it 7) x 10^"*^'^ y^^^- This possibility is 
envisaged in [25] whose authors use d/a = —G/G by speaking about a small 
radial drift of —(6 it 13) x 10~^ m yr~^ in an orbit at 1 AU. 
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3 The evolution of the Earth-Sun system 



In this Section I will not consider other effects which may affect the final 
evolution of the Sun-Earth system like the tidal interaction between the 
Earth and the tidal bulges of the giant solar photosphere, and the drag 
friction in the motion through the low chromosphere [1]. For the Earth, by 
assuming the values a = 1.00000011 AU, e = 0.01671022 at the epoch J2000 
(JD 2451545.0) with respect to the mean ecliptic and equinox of J2000 and 
fi/^ = -9 X 10-1'' yr-\ eq. ([Ml) yields 

Ar(27r) = 1.3 x 10^^ m. (45) 

This means that at every revolution the position of the Earth is shifted 
along the true line of the apsides (which coincides with the osculating one 
because of the absence of perihelion precession) by 1.3 cm. This result is 
confirmed by my numerical integrations and the discussion of Section [21 
indeed, it can be directly inferred from Figure [2] by multiplying the value 
of Ar at t = 1 yr by 9 X 10-^^. By assuming that the Sun will continue 
to lose mass at the same rate for other 7.58 Gyr, when it will reach the tip 
of the RGB in the HR diagram [l], the Earth will be only 6.7 x 10"^ AU 
more distant than now from the Sun at the perihelion. Note that the value 
9 X 10~^'^ yr-i is an upper bound on the magnitude of the Sun's mass loss 
rate; it might be also smaller [Ij like, e.g., 7 x 10" yi"-^ which would yield 
an increment of 5.5 x 10^^ AU. Concerning the effect of the other planets 
during such a long-lasting phase, a detailed calculation of their impact is 
beyond the scope of the present paper. By the way, I wish to note that 
the dependence of Ar(27r) on the eccentricity is rather weak; indeed, it 
turns out that, according to eq. ([26]) . the shift of the perihelion position 
after one orbit varies in the range 1.3 — 1.1 cm for < e < 0.1. Should 
the interaction with the other planets increase notably the eccentricity, the 
expansion of the orbit would be even smaller; indeed, for higher values of 
e like, e.g., e = 0.8 it reduces to about 3 mm. By the way, it seems that 
the eccentricity of the Earth can get as large as just 0.02 — 0.1 [Ml [27l [28] 
over timescales of ~ 5 Gyr due to the N— body interactions with the other 
planets. In Table [T] I quote the expansion of the orbits of the other planets 
of the solar system as well. It is interesting to note that Mercur}0 and 
likely Venus are fated at the beginning of the RGB; indeed, from Figure 2 
of |lj it turns out that the Sun's photosphere will reach about 0.5 — 0.6 AU, 

^It might also escape from the solar system or collide with Venus over 3.5 Gyr from 
now [261 [271 [28]. 



15 



Table 1: Expansion of the orbits, in AU, of the eight planets of the solar 
system in the next 7.58 Gyr for M /M = —9 x 10~^^ yr"^- I have neglected 



Planet 


Ar (AU) 


Mercury 


2 X lO""* 


Venus 


5 X 10"'' 


Earth 


7 X 10-4 


Mars 


9 X lO"'' 


Jupiter 


3 X 10-3 


Saturn 


6 X 10-3 


Uranus 


1 X 10-2 


Neptune 


2 X 10-2 



while the first two planets of the solar system will basically remain at 0.38 
AU and 0.72 AU, respectively, being the expansion of their orbits negligible 
according to Table [TJ After entering the RG phase things will dramatically 
change because in only ~ 1 Myr the Sun will reach the tip of the RGB phase 
loosing mass at a rate of about —2 x 10"^ y^-^ and expanding up to 1.20 
AU In the meantime, according to my perturbative calculations, the 
perihelion distance of the Earth will increase by 0.25 AU. I have used as 
initial conditions for /x, a and e their final values of the preceding phase 7.58 
Gyr- long. In Table [2] I quote the expansion experienced by the other planets 
as well; it is interesting to note that the outer planets of the solar system 
will undergo a considerable increase in the size of their orbits, up to 7.5 AU 
for Neptune, contrary to the conclusions of the numerical computations in 
|29j who included the mass loss as well. I have used as initial conditions the 
final ones of the previous MS phase. Such an assumption seems reasonable 
for the giant planets since their eccentricities should be left substantially 
unchanged by the mutual N-body interactions during the next 5 Gyr and 
more [261 [27t I28j : concerning the Earth, should its eccentricity become as 
large as 0.1 due to the N-body perturbations [261 1271 128] , after about 1 Myr 
its radial shift would be smaller amounting to 0.22 AU. Mutual N-body 
interactions have not been considered, thus hardly preventing our planet 
to escape from engulfment in the expanding solar photosphere. Concerning 
the result for the Earth, it must be pointed out that it remains substantially 
unchanged if I repeat the calculation by assuming a circularized orbit during 
the entire RGB phase. Indeed, it is possible to show that by adopting as 
initial values of o and the final ones of the previous phase I get that 
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Table 2: Expansion of the orbits, in AU, of the eight planets of the solar 
system in the first 1 Myr of the RGB for M/M = -2 x lO"'^ yr'^ I have 
neglected mutual N-body interactions and other phenomena like the effects 
of tidal bulges and chromospheric drag for the inner planets. 



Planet 


Ar (AU) 


Mercury 


7 X 10"^ 


Venus 


1.8 X 10-1 


Earth 


2.5 X 10-1 


Mars 


3.4 X 10-1 


Jupiter 


1.24 


Saturn 


2.25 


Uranus 


4.57 


Neptune 


7.46 



after ~ 1.5 Myr Ar has changed by 0.30 AU. Note that my results are in 
contrast with those in [1] whose authors obtain more comfortable values for 
the expansion of the Earth's orbit, assumed circular and not influenced by 
tidal and frictional effects, ranging from 1.37 AU (|/t/ju| = 7 x yr"^) 
to 1.50 AU = 8 X 10-1"^ yr-i) and 1.63 AU (lA^I = 9 x 10-^^ yr-i). 

However, it must be noted that such a conclusion relies upon a perturbative 
treatment of eq. ([2]) and by assuming that the mass loss rate is constant 
throughout the RGB until its tip; in fact, during such a Myr the term 
{ft/ fj,)At would get as large as 2 x 10-^. In fact, by inspecting Figure 4 of [l] 
it appears that in the last Myr of the RGB a moderate variation of M/M 
occurs giving rise to an acceleration of the order of M/M ~ lO^i^ yr^^- 
Thus, a further quadratic term of the form 

should be accounted for in the expansion of eq. ([2]). A perturbative treat- 
ment yields adequate results for such a phase 1 Myr long since over this 
time span eq. ()46p would amount to ~ 5 x lO-'^. However, there is no need 
for detailed calculations: indeed, it can be easily noted that the radial shift 
after one revolution is ^ 

Ar(2^) oc —. (47) 

After about 1 Myr eq. (j47p yields a variation of the order of 10~^ AU, which 
is clearly negligible. 
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4 Conclusions 



I started in the framework of the two-body Newtonian dynamics by using a 
radial perturbing acceleration linear in time and straightforwardly treated it 
with the standard Gaussian scheme. I found that the osculating semimajor 
axis o, the eccentricity e and the mean anomaly M secularly decrease while 
the argument of pericentre a; remains unchanged; also the longitude of the 
ascending node and the inclination i are not affected. The radial distance 
from the central body, taken on the fixed line of the apsides, experiences a 
secular increase Ar. For the Earth, such an effect amounts to about 1.3 cm 
yr~^. By numerically integrating the equations of motion in Cartesian co- 
ordinates I found that the real orbital path expands after every revolution, 
the line of the apsides does not change and the apsidal period is larger than 
the unperturbed Keplerian one. I have also clarified that such results are 
not in contrast with those analytically obtained for the Keplerian orbital 
elements which, indeed, refer to the osculating ellipses approximating the 
true trajectory at each instant. I applied my results to the evolution of the 
Sun-Earth system in the distant future with particular care to the phase in 
which the Sun, moved to the RGB of the HR, will expand up to 1.20 AU 
in order to see if the Earth will avoid to be engulfed by the expanded solar 
photosphere. My answer is negative because, even considering a small accel- 
eration in the process of the solar mass- loss, it turns out that at the end of 
such a dramatic phase lasting about 1 Myr the perihelion distance will have 
increased by only Ar w 0.22—0.25 AU, contrary to the estimates in [1] whose 
authors argue an increment of about 0.37 — 0.63 AU. In the case of a circular 
orbit, the osculating semimajor axis remains unchanged, as confirmed by a 
numerical integration of the equations of motion which also shows that the 
true orbital period increases and is larger than the unperturbed Keplerian 
one which remains fixed. Concerning the other planets, while Mercury will 
be completely engulfed already at the end of the MS, Venus might survive; 
however, it should not escape from its fate in the initial phase of the RGB 
in which the outer planets will experience increases in the size of their orbits 
of the order of 1.2 - 7.5 AU. 

As a suggestion to other researchers, it would be very important to 
complement my analytical two-body calculation by performing simultaneous 
long-term numerical integrations of the equations of motion of all the major 
bodies of the solar system by including a mass-loss term in the dynamical 
force models as well to see if the N-body interactions in presence of such 
an effect may substantially change the picture outlined here. It would be 
important especially in the RGB phase in which the inner regions of the 
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solar system should dramatically change. 
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